Perfectly Conducting Channel and Universality Crossover in Disordered 

Nano-Graphene Ribbons 
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The band structure of graphene ribbons with zigzag edges have two valleys well separated in mo- 
mentum space, related to the two Dirac points of the graphene spectrum. The propagating modes 
in each valley contain a single chiral mode originating from a partially flat band at band center. 
This feature gives rise to a perfectly conducting channel in the disordered system, if the impu- 
rity scattering does not connect the two valleys, i.e. for long-range impurity potentials. Ribbons 
with short-range impurity potentials, however, through inter-valley scattering display ordinary lo- 
calization behavior. The two regimes belong to different universality classes: unitary for long-range 
impurities and orthogonal for short-range impurities. 

PACS numbers: 72.10.-d,72.15.Rn,73.20.At,73.20.Fz,73.23.-b 



The recent fabrication of graphene devices, combined 
with observation of half- integer quantum Hall effect [l[ 
and the intrinsic 7r-phase shift of the Shubnikov-de Haas 
oscillations 0, has once more ignited an intense dis- 
cussion on this old fascinating system. Due to the 
two-dimensional honeycomb structure, the itinerant tt- 
electrons near the Fermi energy behave as massless Dirac 
fermion. The valence and conduction bands touch coni- 
cally at two nonequivalent Dirac points, called and 
K- points, which possess opposite chirality 

In graphene, the presence of edges can have strong im- 
plications for the spectrum of the 7r-electrons Q • There 
are two basic shapes of edges, armchair and zigzag which 
determine the properties of graphene ribbons. In ribbons 
with zigzag edges, localized states appear at the edge 
with energies close to the Fermi level 01 ■ In contrast, 
edge states are absent for ribbons with armchair edges. 
Recent experiments support the evidence of edge local- 
ized states 0]. The electronic transport through zigzag 
ribbons shows a number of intriguing phenomena such 
as zero-conductance Fano resonances [6| , vacancy config- 
uration dependent transport Q, valley filtering 3 ^-^d 
half- metallic conduction [9|. 

In this Letter, we focus on disorder effects of the elec- 
tronic transport properties of graphene zigzag ribbons. 
The edge states play an important role here, since they 
appear as special modes with partially flat bands and 
lead under certain conditions to chiral modes. There is 
one such mode of opposite orientation in each of the two 
valleys, which are well separated in fc-space. The key re- 
sult of this study is that for disorder without inter-valley 
scattering a single perfectly conducting channel emerges 
associated with such a chiral mode. This mode disap- 
pears as soon as inter-valley scattering is possible. This 
distinction depends on the range of the impurity poten- 
tials. We will show that as a function of the impurity 
potential range a crossover from the orthogonal to the 
unitary universality class occurs which is connected with 
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FIG. 1: (a) Structure of graphene zigzag ribbon. The dis- 
ordered region with randomly distributed impurites lies in 
the shaded region and has the length L. The lattice con- 
stant is a and the ribbon width A'' is defined as the number 
of the zigzag chains, (b) Energy dispersion of zigzag rib- 
bon with N — 10. The valleys in the energy dispersion near 
k — 27r/3a {k — —2n/3a) originate from the Dirac K-^-{K-)- 
point of graphene. The red-filled (blue-unfilled) circles denote 
the right (left)-moving open channel at the energy Eq. In the 
left (right) valley separately the degeneracy between right and 
left moving channels is missing due to one excess right(left)- 
going mode. The time-reversal symmetry under the intra- 
valley scattering is also broken. 



the presence or absence of time reversal symmetry (TRS). 

We describe the electronic states of nanographites by 
the tight-binding model 



H = J2^.,\^){J\+J2V,\^){^\ 



(1) 



where = — 1 if i and j are nearest neighbors, and 



2 



otherwise. \i) represents the state of the p^-orbital on site 
i neglecting the spin degrees of freedom. In the follow- 
ing we will also apply magnetic fields perpendicular to 
the graphite plane which are incorporated via the Peierls 

phase: 7ij ^i^j exp i2TT{e/ch) dl • A , where A is 
the vector potential. The second term in Eq. ([1]) repre- 
sents the impurity potential, Vi = V{ri) is the impurity 
potential at a position rt. 



(b) 




FIG. 2: L-dependence of the average of dimensionless con- 
ductance {g) for zigzag ribbon with A'^ = 10, d/a — 1.5 (no 



inter- valley scattering) , 



1.0, and rii'. 



0.1. 6000 sam- 



ples with different impurity configuration are included in the 
ensemble average. 

As shown in FigUJa), our zigzag ribbons are character- 
ized by the width TV, the number of zigzag chains, and 
L denotes the length of the disordered region. In Fig. 
[TJb), we display the band structure for the zigzag rib- 
bon with A'' = 10. Note that zigzag ribbons are metallic 
for all widths at finite doping because of the presence 
of a partial flat band at zero energy induced by edge 
states. These edge states lead in the clean limit to the 
characteristic conductance odd-number quantization, i.e. 
g = 2n + I as the dimensionless conductance per spin 
(n = 0, ±1,±2,...) d, There are two valleys, at 

k± — ±27r/3, each of which possesses one excess mode 
which violates the balance between the number left- and 
right-moving modes (Figll]). 

In our model we assume that the impurities are ran- 
domly distributed with a density riimp, and the potential 
has a Gaussian form of range d 



V(n) = ^ uexp 

ro {ra7idom) 



In - t-qI 

d2 



(2) 



where the strength u is uniformly distributed within the 
range |u| < um. Here um satisfies the normalization 



condition: um J2ri 



{full space) 



exp I 



}/d') /(V3/2) = uo. 



The range of the impurity potential is crucial for the 
transport properties. Since the momentum difference be- 
tween two valleys is rather large, Ak = k+ — k- = An /3a, 
only short-range impurities (SRI) with a range smaller 
than the lattice constant causes inter-valley scattering. 
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FIG. 3: (a) Distribution of the transmission eigenvalues A: 
p(A), at = 0.5 for L/a = 50,250,1000, with d/a = 2.0. 
-E = 0.5 leads to 3 incident channels. 12,000 samples with 
difi'erent impurity configurations are included in the distri- 
bution, (b) Distribution of the dimensionless conductance g: 
X(A), at L/a = 1000 for the same parameter set. 



Long-range impurties (LRI), in contrast, restrict the 
scattering processes to intra-valley scattering [isj . 

We briefly discuss here the relation between valleys in 
the zigzag ribbons and graphene. The electronic states 
near the Dirac point can be described by the k • p Hamil- 
tonian 



Hk.p = 7 



(3) 



acting on the 4-component pseudo-spinor Bloch functions 
•5 = [4'K+A,(l}K+B,4>K_A,(t>K_B], which characterize 
the wave functions on the two crystalline sublattices (A 
and B) for the two Dirac points (valleys) K±. Here, 7 is 
the band parameter, k^(ky) are wavenumber operators, 
and r° is the 2x2 identity matrix. Pauli matrices a^-y-^ 
act on the sublattice space {A, B), while r^-y-^ on the 
valley space {K±). Since the outermost sites along 
(A^***) zigzag chain are B(A)-sublatticc, an imbalance be- 
tween two sublattices occurs at the zigzag edges leading 
to the boundary conditions 



(4) 



where r[j] stands for the coordinate at i*^ zigzag chain. It 
can be shown that the valley near k = Sir /2a in Fig. 1(b) 
originates from the iir_|_ -point, the other valley at fc = 
-37r/2a from K.-point Q. 

Ignoring the spins the graphene systems have TRS 
with respect to the operator T represented by the com- 
plex conjugation C. Pairs of time reversed states are 
formed across the two- valleys (Dirac points) as expressed 
in the above pseudo-spinor formulation by T = — ?((t^ (g) 
t'^)C, where the A-B-sublattices act as pseudospin de- 
grees of freedom. The boundary conditions which treat 
the two sublattices asymmetrically leading to edge states 
give rise to a single special mode in each valley. Con- 
sidering now one of the two valleys separately, say the 
one around fc = fc_|_ we see that the TRS is effectively 
violated in the sense that we find one more left-moving 
than right-moving mode. If we restrict ourselves to dis- 
order promoting only intra-valley scattering, transport 
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FIG. 4: L-dependence of the averaged dimensionless conduc- 
tance for zigzag ribbons with A'^ = 10, short-ranged impurity 
potentional [d/a = 0.05, inter-valley scattering), uq — 1.0, 
and riimp. ~ 0.1. 5000 samples with different impurity con- 
figurations are in the ensemble average. 



properties would resemble those of a system with a chi- 
ral mode which is oppositely oriented for the two valleys. 
In this sense such a system violates TRS. On the other 
hand, disorder yielding inter-valley scattering would re- 
store this TRS as both valleys together incorporate a 
complete set of pairs of time-reversed modes. Thus we 
expect to see qualitative differences in the properties if 
the range of the impurity potentials is changed. 

In order to demonstrate this we now turn to the dis- 
cussion of the transport properties. The dimensionless 
electrical conductance is calculated using the Landauer- 
Biittiker formula[i3|, giE) = Tt{W), where t{E) is the 
transmission matrix through the disordered region. This 
transmission matrix can be calculated by means of the 
recursive Green function method @, [llj . We focus first 
on the case of LRI using a potential with d/a — 1.5 which 
is already sufficient to avoid inter-valley scattering. Fig[2] 
shows the averaged dimensionless conductance as a func- 
tion of L for different incident energies, averaging over an 
ensemble of 6000 samples with different impurity config- 
urations for ribbons of the width = 10. The potential 
strength and impurity density are chosen to be wq = 1.0 
and riimp. =0.1, respectively. As a typical localization ef- 
fect we observe that (g) gradually decreases with growing 
length L (Figl2]). Surprisingly, the ribbons remain highly 
conductive even at the length oi L = 1500a, i.e. more 
than 350nm in the real system. Actually, (g) converges 
to (g) = 1, indicating the presence of a single perfectly 
conducting channel. It can be seen that {g){L) has an 
exponential behavior as (g) — 1 ^ exp(— L/^) with ^ as 
the localization length^]. 

We performed a number of tests to confirm the pres- 
ence of this perfectly conducting channel. First of all, 
it exists up to L = 3000a for various ribbon widths up 
to iV = 40 for the potential range {d/a = 1.5). More- 
over it remains for LRI with d/a = 2.0,4.0,6.0,8.0, and 
uq = 1.0, nimp. = 0.1 and N = 10. As the effect is con- 
nected with the subtle feature of an excess mode in band 
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FIG. 5: (a) Magnetic field dependence of the inverse local- 
ization length for various incident energies. The fiUed 
symbols indicate data sets of systems with SRI being rather 
sensitive to magnetic fields, and empty symbols denote data 
sets for LRI which are almost insensitive to the field. The 
magnetic flux passing through a single hexagon ring is mea- 
sured in units of the quantum flux {(j}o = ch/e). (b) Enlarged 
view of (a) of LRI data 5000 samples with different configu- 
rations are included in the averages. 

structure, it is natural that the result can only be valid 
for sufficiently weak potentials. For potential strengths 
comparable to the energy scale of the band structure, 
e.g. the energy difference between the transverse modes, 
the result should be qualitatively altered[3|. Deviations 
from the limit (g) 1 also occur, if the incident energy 
lies at a value close to the change between g = 2n — 1 
and g — 2n + 1 for the ribbon without disorder. This is 
for example visible in above calculations for E = 0.4 
where the limiting value (g) < 1 (Fig [2]). As a fur- 
ther test we evaluate the distribution of the transmis- 
sion eigenvalues and dimensionless conductance for fixed 
wire length. In Fig[3]^a), the distribution of the eigenval- 
ues A of the Hermite matrix, tt^ , is depicted for various 
wire lengths. With growing length L a progressive sep- 
aration of the transmission eigenvalues emerges with a 
strong peak close to (localization) and at 1 (perfect con- 
duction channel). The distribution of the conductance g 
(trace of the transmission matrix Tr(ftt))^ jg depicted in 
Figl3fb) for samples in the long- wire limit. Obviously, g 
only distributes above g = 1 with a singularity at 1. 

Turning to the case of SRI the inter-valley scatter- 
ing becomes sizable enough to ensure TRS, such that 
the perfect transport supported by the effective chiral 
mode in a single valley ceases to exist. SRI causes true 
back scattering. For a comparison, we show the ribbon 
length dependence of the averaged conductance in FigH) 
For any incident energy the electrons tend to be local- 
ized and the averaged conductance decays exponentially, 
(g) ~ exp(— ^/L), without developing a perfect conduc- 
tion channel. 

In order to demonstrate that the qualitative difference 
between the two regimes, LRI and SRI, is indeed con- 
nected with TRS, we study the effect of magnetic field 
coupling to the electrons through the Peierls phase. For 
the time reversal symmetric situation resulting from SRI 
scattering the magnetic field removing TRS should have 



4 



a stronger effect than for the case of LRI where TRS 
is broken aheady at the outset. We use the locahza- 
tion length ^ as an indicator. In FigO the field depen- 
dence of the inverse localization length is shown for var- 
ious incident energies (filled symbols for SRI and empty 
symbols for LRI). Indeed the locahzation length displays 
a stronger field dependence than the LRI. Actually for 
LRI even a so-called anti-localization behavior with in- 
creasing field is visible consistent with recent reports on 
graphene Til [itI. . Note that for E < 0.4 only a single 
channel is involved in the conductance such that for LRI 
no localization occurs, i.e. = 0. 

The presence of one perfectly conducting channel in 
disordered quantum wires with symplectic symmetry and 
an odd number of channels has recently been analyzed 
using random matrix theory [l9| . The symplectic sym- 
metry of such systems is based on the skew-symmetry of 
the reflection matrix, *r = — r [23|. A realization can be 
found in the disordered metallic carbon nanotubes with 
LRI. On the other hand, zigzag ribbons without inter- 
valley scattering are not in the symplectic class, since 
they break TRS in the special way. The decisive feature 
for a perfectly conducting channel is the presence of one 
excess mode in each valley. Note that this is in contrast 
to graphene for which each mode has a partner mode 
of reversed velocity in the same valley. For single val- 
ley transport the reflection matrix has a non-square form 



{N^;'''xN^;'' with N^;'''' = nP + 1, where ivf^ (ivi*^ ) is the 
number of the reflection (incident) channels). Recently 
Hirosc et. al. pointed out that non-square reflection 
matrices with unitary symmetry give rise to a perfectly 
conducting channel 21 1 . 

Eventuelly we can identify the universality classes of 
zigzag ribbons. For LRI they belong to the unitary class 
(no TRS), while for SRI with inter- valley scattering they 
are in the orthogonal class (with overall TRS). This clas- 
sification is compatible with the behavior in a magnetic 
field. 

Analogous symmetry considerations can be applied to 
armchair ribbons. In this case the two valleys merge into 
a single one at A; = 0. TRS is conserved irrespective of 
the impurity potential range, if there is no magnetic field. 
Consequently, disordered armchair ribbons belong always 
to the orthogonal class and do not provide a perfectly 
conducting channel. In view of the fact that graphene is 
known to be symplectic (orthogonal) for LRI (SRI) [13], 
it is quite intriguing to realize that the edges influence the 
universality class, as long as the phase coherence length 
is larger than the system size of nanographenes. 

The unusual energy dispersion due to their edge states 
gives rise to the unique property of zigzag ribbons. Con- 
cerning transport properties for disordered systems the 
most important consequence is the presence of a per- 
fectly conducting channel. The origin of this effect lies in 
the single-valley transport which is dominated by a chi- 
ral mode. On the other hand, large momentum transfer 



through impurities with short-range potentials involves 
both valleys, destroying this effect and leading to usual 
Anderson localization. The obvious relation of the chiral 
mode with time reversal symmetry leads to the classifi- 
cation into the unitary and orthogonal class depending 
on the range of impurity potential. Since the inter- valley 
scattering is weak in the experiments of graphene, we 
may assume that these conditions may be realized also 
for ribbons. Naturally defects in the ribbon edges and 
vacancies would be rather harmful for the experiment 
making this type of experiment very challenging 
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